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Method for designing a nuclear fuel assembly having a 
damping guide tube, a corresponding system, computer 

programme and medium 

The present invention relates to a method for designing a 
nuclear fuel assembly which is intended to be positioned 
in a nuclear reactor, the assembly comprising a plurality 
of guide tubes, and a control cluster which itself 
comprises a plurality of control rods and a support for 
control rods, the control rods and the guide tubes 
extending in parallel with a longitudinal direction, each 
control rod being received in a guide tube in order to 
form pairs comprising guide tubes/control rods, each 
guide tube comprising a lower damping portion which 
comprises at least a portion of reduced inside diameter, 
which portion is intended to contain a fluid for damping 
the fall of the control rod which is received in the 
guide tube, the portion of reduced inside diameter 
.surrounding the control rod with a radial passage gap 
when the control rod enters the guide tube. 

It will be appreciated that nuclear fuel assemblies must 
be dependable in order to allow reliable operation of 
nuclear reactors. 



Thus, design and construction provisions for such 
assemblies have been drawn up. 

These provisions impose a general framework and minimum 
criteria which the assembly constructors must take into 
consideration. 

As far as the guide tubes are concerned, the design 
provisions require verification that they remain stable 
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under axial compression and that the mechanical integrity 
thereof is not affected during such compression. These 
criteria are aimed in particular at taking into 
consideration the axial compression of the guide tubes 
5 which results from contact against the upper bearing 
plate of the core by springs carried by an upper end 
piece of the assembly. 

Although the criteria imposed by the design provisions 
10 allow assemblies to be designed with satisfactory 
reliability, these are merely minimum criteria and it is 
desirable to develop a method which allows even more 
reliable assemblies to be designed. 

15 For all that, it is also desirable to reduce the safety 
margins during design in order to reduce the mass and the 
cost of the assemblies constructed. 

An object of the invention is to overcome this problem by 
20 providing a method which allows nuclear fuel assemblies 
which are more reliable to be designed, whilst limiting 
the design margins. 

To this end, the invention relates to a method for - 
25 designing a nuclear fuel assembly which is intended to be 
positioned in a nuclear reactor, the assembly comprising 
a plurality of guide tubes, and a control cluster which 
itself comprises a plurality of control rods and a 
support for control rods, the control rods and the guide 
30 tubes extending in parallel with a longitudinal 
direction, each control rod being received in a guide 
tube in order to form pairs comprising guide 
tubes/control rods, each guide tube comprising a lower 
damping portion which comprises at least a portion of 



reduced inside diameter, which portion is intended to 
contain a fluid for damping the fall of the control rod 
which is received in the guide tube, the portion of 
reduced inside diameter surrounding the control rod with 
5 a radial passage gap when the control rod is introduced 
in the guide tube, characterised in that the method 
comprises, for at least one pair comprising a guide 
tube/control rod, the following steps: 

10 a) establishing the falling speed of the control rod 
upon entry into the lower damping portion when the 
control cluster falls in the event of a shutdown of the 
nuclear reactor, 

b) establishing, based on the falling speed established 
15 in step a) , the progression of the falling speed of the 

control rod in the lower damping portion, 

c) establishing, based on the progression of the speed 
established in step b) , a maximum elevated pressure 
produced in the liquid contained in the lower damping 

20 portion, and 

d) establishing, based on the maximum elevated pressure 
established in step c) , a maximum circumferential stress 
produced in the lower damping portion. 

25 According to specific embodiments, the method can 
comprise one or more of the following features, taken in 
isolation or according to all technically feasible 
combinations : 

30 - the method further comprises a step for verifying, 
using the maximum circumferential stress established in 
step d) , that a maximum stress admissible by the guide 
tube has not been exceeded. 
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the establishing step b) is carried out using a higher 
value for the radial passage gap and the establishing 
step c) is carried out using a lower value for the radial 
passage gap, 

5 - the higher value is a maximum statistical value for the 
passage gap, 

the lower value is a minimum statistical value for the 
passage gap, 

the support of the control cluster comprising a helical 
10 spring for damping the impact of the support against an 
upper end piece of the assembly in the event of the 
control cluster falling during a shutdown of the nuclear 
reactor, the method further comprises the following 
steps : 

15 e) establishing the progression of the speed of the 

control cluster after the impact of the support against 
the upper end piece, 

f) establishing, based on the progression of the speed 
established in step e) , a maximum longitudinal load for 

20 compression of the spring, and 

g) establishing, based on the maximum longitudinal load 
for compression, at least a maximum shearing stress in 
the spring, 

a maximum shearing stress is a shearing stress along 
25 the neutral axis of the spring, 

a maximum shearing stress is a shearing stress along 
the axis of the spring nearest the longitudinal centre 
axis thereof, 

the method further comprises a step for verifying, 
30 using a maximum shearing stress established in step g) , 
that a maximum stress admissible by the spring has not 
been exceeded. 
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The invention further relates to a system for designing a 
nuclear fuel assembly, characterised in that it comprises 
means for carrying out the steps of a method as defined 
above . 

5 

According to a variant, the system comprises a computer 
and storage means, in which at least a programme 
comprising instructions for carrying out steps of the 
method for designing a nuclear fuel assembly is stored. 

10 

The invention further relates to a computer programme 
comprising instructions for carrying out the steps of a 
method as defined above. 

15 The invention also relates to a medium which can be used 
in a computer and on which a programme as defined above 
is recorded. 

The invention will be better understood from a reading of 
20 the description below which is given purely by way of 
example with reference to the appended drawings, in 
which : 

Figure 1 is a schematic, perspective cut-away view of a 
25 nuclear fuel assembly which is designed by a method 
according to the invention. 

Figure 2 is a schematic, partially sectioned side view 
drawn to an enlarged scale of the structure of the spider 
of the assembly of Figure 1, 
30 - Figure 3 is a partial schematic side view of the 
assembly of Figure 1, illustrating more particularly a 
pair comprising a guide tube/control rod. 

Figure 4 is a block diagram illustrating the system for 
designing the assembly of Figure 1, 
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Figure 5 is a flow chart illustrating successive steps . 
of the design method carried out by the system of Figure 
4, 

Figure 6 is a progression curve of the falling speed of 
5 a control rod before it is introduced in the lower 
portion of the corresponding guide tube, this progression 
being calculated by the system of Figure 4, and 

Figure 7 is a progression curve of the falling speed of 
the same control rod in the lower portion of the 
10 corresponding guide tube, this progression being 
calculated by the system of Figure 4 . 

Figure 1 illustrates a nuclear fuel assembly 1 which 
mainly comprises a square-based lattice 2 for nuclear 
15 fuel rods 3 and a control cluster 4 . 

The assembly 1 comprises grids 5 for maintaining the rods 
3, which grids 5 are distributed over the height of the 
rods 3 • A lower end piece 6 is arranged under the lower 
20 ends of the rods 3 and an upper end piece 7 above the 
upper ends of the rods 3 . The upper end piece 7 is 
provided with springs 8 for pressing against the upper 
bearing plate of the reactor core, in which the assembly 
1 is intended to be placed . 

25 

The control cluster 4 comprises a plurality of control 
rods 10, for example, 24. Conventionally, the control 
rods 10 comprise a material which absorbs neutrons. 

30 The rods 3 and 10 extend in parallel with a vertical 
longitudinal direction L. 



The rods 10 are carried at the upper ends thereof by a 
support 11 which is generally referred to as a spider. 
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As illustrated more particularly in Figure 2, the spider 
11 comprises a vertical central upper head 12 and a 
series of arms or vanes 13 which extend radially outwards 
5 from the lower end of the upper head 12 as far as the 
radially outer ends 14 thereof. 

Each control rod 10 is connected to an arm 13 at the 
upper end thereof • 

10 

The upper head 12 of the spider 11 has a central blind 
hole 15 which opens towards the bottom and in which a 
damping helical spring 16 is received. The spring 16 
extends vertically along a centre axis A. A tightening 
15 screw 17 extends substantially over the entire height of 
the hole 15 and is screwed into the wall 18 delimiting 
the upper portion of the hole 15. 

The lower portion of the screw 17 extends through the 
20 base of a retaining ring 2 0 which rests on the lower end 
of the spring 16. The head 21 of the screw 17 rests, at 
the top, against the base of the retaining ring 2 0 in 
order to press the spring 16 against the wall 18 of the 
upper head 12 . 

25 

As illustrated in Figure 3 for a control rod 10, each 
control rod 10 is received in a respective guide tube 24 
which is arranged in the lattice 2 of fuel rods 3 . In 
this manner, 24 pairs comprising a guide tube/control rod 
30 are formed. Since each of these pairs has a similar 
structure, only one will be described below. 

The guide tube 24 extends from the lower end piece 6 as 
far as the upper end piece 7. The guide tube 24 comprises 
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a lower portion 2 6 of reduced inside diameter and an 
upper portion 27. The lower portion 26 is connected to 
the lower end piece 6 by a collared screw 28, through 
which a vertical through-hole 29 extends. 

5 

The lower portion 2 6 of the guide tube 24 surrounds the 
control rod 10 with a radial passage gap J. 

The upper portion 27 is fixed to the upper end piece 7 
10 and opens at the outside of the assembly 1. 

Lateral apertures 30, only one of which can be seen in 
Figure 4, are provided in the upper portion 27 near the 
lower portion 26. 

15 

When the assembly 1 is placed in a nuclear reactor, the 
cooling liquid of the reactor fills the interior of the 
guide tube 24 . 

20 Conventionally, the control cluster 4 can be moved 
vertically relative to the remainder of the assembly 1 in 
order to allow adjustment of the reactivity during normal 
operation of the reactor, and therefore variations in 
power from zero power up to maximum output depending on 

25 the vertical introduction of the control rods 10 in the 
lattice 2 of rods 3 . The vertical displacement of the 
control cluster 24 is conventionally carried out by way 
of a drive rod which is connected to the upper end of the 
upper head 12 . 

30 

When the reactor is shut down, the drive rod and the 
assembly 4 fall owing to gravity. 
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At the start of this falling movement, the control rods 
10 are guided only by the upper portions 2 7 of the guide 
tubes 24 and have not yet reached the lower portions 26. 

5 Once the falling action has ended, the lower ends of the 
control rods 10 are introduced in the lower portions 26. 
The cooling fluid contained in the portions 2 6 is then 
violently forced, on the one hand, upwards thereby and, 
on the other hand, downwards through the apertures 2 9 of 
10 the collared screws 28. 

Each lower portion 2 6 therefore behaves in the manner of 

a hydraulic damper braking the falling movement of the 
corresponding control rod 10, and therefore of the 
15 assembly 4 . 

This braking phase ends at the end of the travel path 
with the impact of the spider 11 against the upper end 
piece 7 of the assembly 1. 

20 

This impact is carried out by means of the retaining ring 
20. During this impact, the spring 16 is compressed 
vertically in order to absorb the shock. 

25 According to the invention, the assembly 1 has been 
designed in order to take into consideration the specific 
stresses brought about in the assembly by the fall of the 
control cluster 4 during such a shutdown of the reactor. 

30 In this manner, in order to design the assembly 1, in 
particular a data-processing system 32 has been used, as 
illustrated schematically in Figure 4. 



10 

This system 32 comprises, for example, a computer or data 
processing unit 34 comprising one or more processors, 
storage means 36, input /output means 38, and optionally 
display means 40. 

5 

Instructions which can be carried out by the computer 34 
are stored in the form of one or more programmes in the 
storage means 36. 

10 These instructions are, for example, instructions in 
FORTRAN programming code . 

These various instructions, when they are carried out by 
the computer 34, allow the method illustrated by the flow 
15 chart of Figure 5 to be carried out . 

In a first step illustrated by the box 42 of this Figure, 
the computer 34 calculates, based on data 43, the 
progression of the falling speed of a control rod 10 in 
20 the upper portion 2 7 of the corresponding guide tube 24 
in the event of a shutdown of the reactor. 

This calculation can be carried out assuming, for 
example, that the control rod 10 is first subjected to 
25 constant loads : 

gravitational force: fg = Mg, 
Archimedes ' thrust": ^f a = -p gV, 
pressure difference in the core: fc, and 
mechanical friction: fm, 
30 where M and V are the mass and the volume, respectively, 
of the assembly 4 and the drive rod thereof. 



The control rod 10 is also subjected to loads as a 
function of the speed or position thereof, for example. 
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hydraulic friction which can be obtained from: fh = -cl 
(M + pV) v^, with V = speed of the assembly 4 and 
therefore of the rod 10 in question. 

5 Thus, the equation of the movement of the rod in the 
upper portion 27 of the guide tube 24 is as follows: 

(M + pV) 4^ = Sf 
at 

10 This gives: 

4^ =c2-cl 
dt 

with cl = hydraulic friction in the guide tube and 

15 Cl and c2 are, for example, experimental data measured 
during drop tests of the control cluster 4 . These data 
are, with the other data necessary for the calculation, 
such as the mass and the volume of the assembly 4 and the 
drive rod thereof, introduced, for example, in the form 

20 of a file 43 by way of the input/output means 38. 

The computer 34 resolves the equation of the movement of 
the control rod 10, for example, using the NEWTON method. 

25 Thus, the progression of the speed of the control rod 10 
in the upper portion 27 is known as a function of time. 
The profile established in this manner can be displayed 
in the form of a curve by the display means 40. This 
curve is illustrated by Figure 6. 

30 

In this manner, at the end of the step illustrated by the 
box 42, the speed of the control rod 10 is known at the 
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point of entry to the lower damping portion 2 6 of the 
guide tube 24 . 

Based on the results of the step of box 42, the computer 
34 calculates the progression of the speed of the control 
rod 10 during its fall in the lower damping portion 26. 

This step is schematically illustrated by the box 44 . 

This step can be carried out using the following 
equation : 



SCA = cross-section of the rod 10 and 

NCA = number of rods 10 in the assembly 4. 

Therefore, the hypothesis that fc and fm are negligible is 
applied here. 

The difference AP represents the elevated pressure 
produced in the cooling liquid contained in the guide 
tube 24, that is to say, the pressure thereof between the 
lower end of the rod 10 and the pressure present in the 
upper portion 27 of the guide tube 24. 

AP can be established by the following formula: 




with 



c2 = 



fg^ fa ^M-pV 
M+pV'M+pV 



g 





with SM = cross-section of the lower portion 26, 
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SACM = SM - SCA = cross- section of the annular space 
between the rod 10 and the lower portion 26, 

SACTG = STG - SCA, where STG is the cross-section of the 
upper portion 27 of the guide tube 24, 

5 CONTRA =0.4 (l-^^ . 

FECR = coefficient of loss of load owing to friction in 
the lower portion 26, 

"^'^ -Wlwi' 

DM = mean diameter of the guide tube 24 in the upper 
10 portion 27, 

z = height of the rod 10 introduced in the lower portion 
26 of the guide tube 24, and 

Q = fraction of liquid flowing upwards out of the lower 
portion 26. 

15 

The resolution of the equations governing the movement of 
the rod 10 after entry into the lower portion 26 is 
carried out by the computer 34, for example, using the 
RUNGE-KUTTA method, 

20 

Thus, at the end of the step 44, the progression of the 
speed of the control rod 10 in the lower portion 26 of 
the guide tube 24 is known before the impact of the 
spider 11 on the upper end piece 7. 

25 

The speed profile established in this manner can be 
displayed, for example, by the means 40, as illustrated 
in Figure 7. On the curve in Figure 7, the speed profile 
established during step 44 is the portion located to the 
30 left of the point 45. 
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The computer 34 then carries out, in the step of box 46, 
the calculation of the maximum elevated pressure produced 

APmax. 

5 This calculation ' can be carried out, for example, based 
on the formula: 

AP = ^ pQ'v'iEXPA^CONTRA+FECR^CISAxz) . 

The computer 34 carries out, in the step of box 48, the 
10 calculation of a circumferential stress and maximum 
normal aeMAX/ to which the lower portion 2 6 of the guide 
tube 24 is subjected owing to the maximum elevated 
pressure APmax. 

15 This stress can be calculated based on the formula: 

where DPM = inside diameter of the lower portion 2 6 and 
20 EMP = minimum thickness of the wall of the lower portion 
26 . 

The system 32 can then provide, owing to the input /output 
means 38, a first result in the form of a file 49 
25 containing the value aoMAx established, and optionally the 
maximum elevated pressure APmax established, 

Next, the system 32 carries out the calculation of the 
progression of the speed of the control rod 10 after it 
30 comes into contact with the spider 11 and the upper end 
piece 7 . 



is 

This calculation step is illustrated by the box 50 in 
Figure 5 . 

This calculation can be carried out, for example, using 
5 the following equation when the ring 20, and therefore 
the spider 11, is in contact with the upper end piece 7: 

(M + pV) % = (M - pV) g - PRCH - K (z- LAI) - c3 v 
at 

with PRCH = preload of the spring 16 = PRCMP x K, where 
PRCMP is the precompression of the spring 16 and K the 
10 rigidity of the spring 16, 

LAI = distance travelled by the control rod in the lower 
portion 2 6 before impact, and 

c3 = coefficient of hydraulic damping in order to model 
the damping in the lower portion 26. 

15 

In the event of a rebound, that is to say, when the 
spider 11 is no longer in contact with the upper end 
piece 7, the equation for movement of the control rod 10 
in question is written as follows: 

20 (M + pV) ^ = (M - pV) g - c3 V 

These two equations are integrated by the computer 34, 
for example, using the RUNGE-KUTTA method. 

Therefore, the step 50 allows the kinematics of the 
25 control cluster 4 to be established during the mechanical 
damping of the shock by the spring 16. The speed profile 
established in this manner can be displayed, for example, 
by the means 40. This profile corresponds to the portion 
located to the right of the point 45 on the curve in 
30 Figure 7 . 

Based on the results of this step, the system 32 carries 
out, in the step 52, the calculation of a maximum 
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20 



25 



vertical compression force Fmtoc/ to which the spring 16 is 
subjected during the mechanical damping. 

This calculation can be carried out, for example, based 
on the following formula: 
Fmax = max {K(z-LAI) + PRCH} 

The system 32 then carries out, in the step of box 54, 
the calculation of an approximate maximum shearing stress 
Tmax in the spring 16 : 



Tmax 5 — 

TtDFR 

with DFN = DER-DFR and 

DER = outside diameter of the spring 16, 
15 DFR = diameter of the wire of the spring 16. 

Subsequently, the system 32 can optionally carry out, 
based on the maximum stress Tmax/ the calculation of 
maximum corrected stresses . 



These stresses can be calculated by multiplying Tmax by 
different factors . 

Thus, it is possible to calculate: 
tmaxi = Tmax x Kc, and 

f MAX2 = T|S4AX X K, 



-A 

c 



with Kc = 1 + 
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The stress Tmaxi corresponds to the shearing stress along 
the neutral axis FN (Figure 2) of the spring 16. The 
stress Tmax2 corresponds to the stress along the axis F2 
(Figure 2) of the spring 16 nearest the vertical centre 
axis A of the spring 16 (see Figure 2) . 

At the end of this step illustrated by the box 56, the 
system 32 provides the various maximum shearing stresses 
calculated, for example, in the form of data stored in a 
file 57, which are transmitted by the input/output means 
38. 

Based on the data contained in the files 49 and 57, which 
have also been stored in the storage means 36, the 
computer 34 will verify that the maximum stresses 
calculated are indeed acceptable for the materials which 
respectively constitute the guide tube 24 and the helical 
spring 16. 

This step has been schematically illustrated by the box 
58 in Figure 5. During such a step, the system 32 will, 
for example, verify that the maximum shearing stresses 
calculated during the steps 54 and 56 are less than 
maximum values admissible by the material which 
constitutes the spring 16. This verification is carried 
out by a comparison of Tmax# ^maxi and Tmax2 with a maximum 
value admissible by the material of the spring 16. 

As far as the maximum circumferential stress aeMAx is 
concerned, the verification can be carried out based on a 
formula of the type: 

f (CTeMAx) < GT 

admissible 
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where aadmissibie refers to the material which constitutes 
the lower portions 2 6 of the guide tubes 24. 

The function f can be a function which takes into 
5 consideration other stresses to which the guide tubes 24 
can be subjected. Such a stress can be a vertical 
compression stress Gp,, to which the guide tubes 24 are 
subjected during the contact of the springs 8 of the 
upper end piece 7 against the upper bearing plate of the 
10 core in order to counterbalance the hydrostatic thrust 
during operation. 

Thus, the function f can be, for example, in the form of 

f (CTGMAXf CJa) = aeMAX + C^A • 

15 

It will be appreciated that this last step, illustrated 
by the box 58, can be carried out by separate software 
which generally carries out the validation of various 
design parameters of the assembly 1 based on results 
20 provided by various pieces of software each dedicated to 
taking into consideration specific operating conditions 
and which include the software which carries out the 
steps 42, 44, 46, 48, 50, 52, 54 and 56. 

25 In general terms, the file 43 comprising the data 43 used 
by the method for the various calculations can comprise 
the data of Table 1 below. 
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Table 1 



Similarly, the file 49 comprising the results from step 
4 8 can comprise the data of Table 2 below. 



APmax: maximum elevated pressure in lower portion 26 


(Pa) 


Zmax: corresponding penetration in lower portion 26 


(m) 


cJeMAx: maximum stress in lower portion 26 


(Pa) 


finax: maximum force on lower end piece 6 


(N) 


tdur: duration of fall in lower portion 26 before impact 


(s) 


vfin: speed of impact of assembly 4 on upper end piece 7 


(m/s) 



Table 2 



The file 57 comprising the results of step 56 can itself 
contain the data of Table 3 below. 
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Fmax^ maximum compression force on spring 16 

Hmax- maximum deflection of spring 16 

'^MAx' approximate maximum stress in spring 

tmaxi ' approximate maximum stress corrected by Kc 

XMAX2: approximate maximum stress corrected by K (Wahl 

coefficient) 

Table 3 




It has been possible to verify by experiment that the 
maximum elevated pressures and the maximum stresses 
5 obtained by means of steps 42, 44, 46 and 48 were 
reliable. In this manner, the first corresponding part of 
the method allows reliable guide tubes 24 to be designed. 
Furthermore, this first part calculates only a single 
stress which appears to be the pertinent stress for the 
10 conditions being considered. Consequently, this first 
part of the method allows the security margins to be 
limited during design, and therefore assemblies which are 
relatively light and economical to be designed. 

15 The second part of the method, which corresponds to steps 
50, 52, 54 and 56, also allows maximum stresses to be 
reliably calculated, as confirmed by experiment. 



Thus, the second part of the method allows a reliable 
20 design to be arrived at by calculation for the spider 
springs 16, which design is found to be advantageous in 
comparison with the method of tests alone which is 
currently imposed by provisions. It will be appreciated 
that the second part of the method calculates only the 
25 small number of stresses, and in particular those on the 
axis F2 of the spring 16 nearest the centre axis A of the 
spring, which are found to be pertinent to the conditions 
envisaged. In this manner, the second part of the method 
allows the design margins to be reduced. 
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In more general terms, the steps 42, 44, 4 6 and 48, on 
the one hand, and 50, 52, 54 and 56, on the other, can be 
carried out by separate pieces of software. 

5 In order to increase the reliability of the calculation, 
for carrying out the first part of the method it is 
possible to use, as the passage gap J, the nominal value 
of the gap, or this nominal value corrected by the 
manufacturing tolerance, or a value resulting from 

10 statistical studies of the distribution of passage gaps J 
obtained in constructed assemblies. 

In a variant, it is possible to use a gap value J which 
is greater for steps 42 and 44 and a smaller gap value J 

15 for steps 46 and 48. This allows a high stress value aeMAx 
to be calculated because the speed reached during the 
fall of the rod 10 in question is high and the volume 
available in the lower portion 26 for the liquid during 
damping is small. However, this high stress value is not 

20 unrealistic and therefore does not lead to unjustified 
design margins, as illustrated by the following example. 

According to a specific variant, the upper value can be a 
maximum value for gap J which is verified with a given 

25 probability, for example, 95%, in constructed assemblies, 
and the lower value can be a minimum value obtained with 
the same probability. This variant allows an 
approximation of the situation where a single pair 
comprising a guide tube/control rod has minimum gap J, 

30 where the maximum stress aeMAx would be reached, and where 
all the other pairs comprising a guide tube/control rod 
have the maximum passage gap J, which would be the most 
extreme case. 
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In some variants, the first part of the method could also 
take into consideration forms of the lower damping 
portion 26 which are different from those described 
previously- In this manner, these lower damping portions 
5 could have a plurality of successive portions of reduced 
diameter, optionally separated by portions of increased 
diameter, generally referred to as cavities. In some 
variants, the first part of the method is carried out 
with collared screws 2 8 which are not perforated by holes 
10 29 • 

In still more general terms, the first part and the 
second part of the design method described can be used 
independently of each other. In this manner, it is 
15 possible to carry out the second part relating to the 
design of the spring 16 without referring to the 
calculation of the elevated pressure AP and the stress 



